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Bulk BiMnOs is the only transition-metal perovskite oxide that is insulating and shows strong 
ferromagnetism. This distinctive behavior would make it a promising candidate as a magnetoelec¬ 
tric multiferroic if it was also a polar material, but experiments have shown that bulk BiMnOs has 
either a very small polarization (below 0.1 pC/cm 2 ) or, most likely, that it is a paraelectric. There 
is also experimental evidence that the polarization in BiMnOs films grown on SrTiOs can be as high 
as 20 pC/cm 2 . Despite of the interest of these behaviors, the diagram of BiMnOs as a function of 
epitaxial strain has remained largely unexplored. In this article, we use first-principles to predict 
that both under enough compressive and tensile epitaxial strain BiMnOs films are ferroelectric with 
a giant polarization around 100 pC/cm 2 . The phases displayed by the films are similar to those 
experimentally found for BiFeOs in similar conditions—at compressive strains, the film is superte- 
tragonal with a large component of the polarization pointing out of plane, while at tensile strains 
the polarization points mostly in plane. Like in BiFeOs films, these phases are antiferromagnetic— 
the orbital ordering responsible for ferromagnetism in BiMnOs is absent in the polar phases. Our 
calculations also show that the band gap of some of these BiMnOs films is substantially smaller than 
gaps typically found in ferroelectric oxides, suggesting it may be a suitable material for photovoltaic 
applications. 

PACS numbers: 77.84.-s, 75.85.+t, 71.15.Mb 


I. INTRODUCTION 


BiMnOs is the only strong ferromagnetic insulator 
among the transition-metal perovskite oxides, a family of 
functional materials whose members display many differ¬ 
ent properties of technological interest. The Curie tem¬ 
perature of BiMnC >3 is around 105 below which sev¬ 
eral groups have reported measurements of its magnetic 
moment close to 4 pb per Mn atom^^ The unique ferro¬ 
magnetic behaviour of BiMnOs is related to the presence 
of orbital ordering of the Mn +3 ions (3d 4 )—elongated 
half-filled d 2 orbitals point towards empty d x 2_ y 2 orbitals 
of neighboring Mn cations^ For many years BiMnOs 
was supposed to have also a weak switchable polariza¬ 
tion (below 0.1 pC/ cm 2 ), and the polar space group C 2 
was assigned to its crystal structure££ However, recent 
electron diffraction experiments by Belik and coworkers^ 
point to a C2/c space group, which is not polar; these re¬ 
sults have been confirmed by neutron powder diffraction 
experiments^ Despite of this, when BiMnOs is grown 
as a epitaxial film on SrTiOs it develops a strong po¬ 
larization, which Jeen and coworkers^ measured to be 
23 pC/ cm 2 ; other groups reported similar results—a po¬ 


larization between 9 and 16 pC/ cm 2 p and signatures of 
ferroelectricity^ 

Perovskite oxides display dramatic changes in their 
properties when grown as epitaxial films on a substrate, 
with the misfit strain imposed by the substrate playing 
a major role in these changesAi Given that other non¬ 
polar perovskites develop ferroelectricity under epitaxial 
strain, it is somewhat surprising that the epitaxial phase 
diagram of BiMnOs is largely unexplored in the search for 
possible coexistence of magnetic and polar orderings. Us¬ 
ing first-principles calculations, Hatt and Spaldin 12 con¬ 
cluded that BiMnOs remains non-polar when compres¬ 
sive or tensile strain is applied to the film; they used the 
C2/c bulk phase as the starting point of their calcula¬ 
tions, subjecting the simulation cell to constraints that 
mimic the strained film. No other phases seem to have 
been studied using computational approaches, although 
earlier calculations by Hill and Rabei^ and by Seshadri 
and Hillii pointed to the existence of ferroelectric insta¬ 
bilities in this material. In this study we explore struc¬ 
tures that are good candidates to become the ground 
state of epitaxial BiMnOs films. At high compressive 
strains, we have looked at the supertetragonal phases 
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with giant polarization that are found in BiFeOs hlms^ 
and in BiCoC >3 bulk and films' these phases are called 
supertetragonal, or T, because they have a large c/a ra¬ 
tio (some of them have lower symmetry than tetragonal). 
At high tensile strains, we have explored phases of low 
energy in BiFeOa (including the Pnma structure of both 
bulk BiFeC >3 and bulk BiMnC >3 under pressure^ and epi¬ 
taxial polar phases of BiFeC >3 under tensile strains). We 
have found that BiMnC >3 and BiFeC >3 films display simi¬ 
lar structures under high tensile and compressive strains; 
therefore, we predict that BiMnC >3 has a strong polariza¬ 
tion when grown in these conditions. 


II. METHODS 


Our first-principles calculations are based on density- 
functional theory (DFT)i 18 ' 19 However, BiMnOs con¬ 
tains highly localized d orbitals that are not described 
accurately within pure DFT. One approach towards a 
better description of the electronic properties of these 
solids is to add a “Hubbard U” term to the energy of the 
system that favors localization on those electrons; this 
requires picking a value of U that reproduces some set of 
experimental or more accurate theory results. Another 
approach is to use a hybrid density functional that in¬ 
cludes a portion of exact Hartree-Fock exchange^ this 
in general ameliorates the problems with d or / elec¬ 
tron delocalization, predicts band gaps for solids that are 
much closer to experimental results than those of pure 
DFT,— and performs better than the “Hubbard U” ap¬ 
proach in perovskite oxides such as BiFeC> 3 . 22 Both ap¬ 
proaches are implemented in VASPr^ the first-principles 
code that we have used to carry out the calculations 
presented in this work. To optimize our typical struc¬ 
tures with Vasp the second approach requires a hun¬ 
dred times more computer time than the first one, so 
we have used the code in the following way: for our ex¬ 
ploratory calculations we applied a Hubbard U follow ing 
the rotationally invariant method described in Ref. l24l . 
and for our final results we used a hybrid functional ac¬ 
cording to the HSE06 prescription. 21 ' As in our previous 
article on BiFeOs,— we worked with the Perdew-Burke- 
Ernzerhof DFT exchange-correlation functional adapted 
to solids (PBEsol) — We used the projector augmented- 
wave method to represent the ionic cores,— solving for 
the following electrons: Mn’s 3 p, 3d, and 4s; Bi’s 5 d, 6s, 
and 6p; and O’s 2s and 2 p. We represented wave func¬ 
tions in a plane-wave basis set truncated at 500 eV. We 
performed integrations within the Brillouin zone by using 
fc-point grids with densities similar to that of a 6 x 6 x 6 
grid for a 5-atom perovskite unit cell. 


III. RESULTS AND DISCUSSION 
A. Metastable Phases of Bulk BiMnOa 

We started our search for metastable phases of 
BiMnC >3 by doing PBEsol+f/ ionic relaxations for the 
seven lowest-energy crystal structures found for bulk 
BiFeC >3 in our previous worh2£ (which can accommodate 
to a 2 x 2 x 2 pseudocubic 40-atom cell) and for the ex¬ 
perimental C2/c ground state configuration of BiMnOs 
(using 40-atom cells that are not pseudocubic). For 
each of these eight configurations we prepared four types 
of magnetic arrangments, as in Ref. |25|—ferromagnetic 
(FM), and antiferromagnetic of the A, C, and G types 
(A-AFM, C-AFM, G-AFM, respectively). We then did 
three types of searches for local minima of the energy: 
(1) we directly relaxed the structures until forces and 
stresses were close to zero; (2) we did a few steps of 
molecular dynamics in order to break possible spureous 
symmetries, and then we relaxed the resulting structures 
until forces and stresses were close to zero; and (3) we 
took the lowest-energy magnetic ordering found so far 
for each structure type and relaxed the atoms imposing 
each of the other three magnetic orderings. At the end, 
we chose the lowest-energy phase for each type of struc¬ 
ture found and for each type of magnetic ordering. In all 
the optimization calculations of this work the final forces 
were below 0.015 eV/A and the final stresses are below 
0.0005 eV/A 3 . 

The process just described lead to the identification 
of the configurations whose energies are given in Fig. Q] 
(top). The corresponding phases are labeled in the fol¬ 
lowing way (directions are given in the pseudocubic set¬ 
ting): the ground state with C2/c symmetry is labeled 
as GS; the paraelectric phase with Pnma symmetry is 
labeled as p- the ferroelectric phase derived from the R3c 
phase that is the ground state of BiFeOs is labeled as 
Raac (since Mn +3 is a d 4 Jahn-Teller active ion, the orig¬ 
inal R3c phase distorts into this one, which has a polar¬ 
ization with a component along [110] and another one 
along [001], corresponding to the Cc monoclinic space 
group); the other phases are supertetragonal T phases 
like the ones mentioned earlier. Three of these T phases 
are local minima of the energy according to our analysis 
of the corresponding force-constant matrices: T aac (orig¬ 
inating from the Cc phase, with a small component of 
the polarization along [110] and a large one along [001]), 
T c (originating from the Pna2\ phase, with polarization 
along [001]), and T' aac (originating from the Pc phase, 
with a small component of the polarization along [110] 
and a large one along [001]). We also found that two of 
these T phases are not local minima of the energy: T ac 
(originating from the Cm phase, with a small compo¬ 
nent of the polarization along [010] and a large one along 
[001]) and T' c (the simplest P4mm tetragonal configura¬ 
tion with polarization along [001]). 

The calculations mentioned in the previous paragraph 
were done using U = 4 eV and J = 1 eV, since these val- 
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TABLE I. Properties of the bulk BiMnOa phases correspond¬ 
ing to minima of the energy found in this work: name given 
to the phase (and most favorable magnetic ordering), space 
group, energy (in meV per formula unit with respect to the 
ground state), band gap (in eV), c/a ratio, and polarization 
vector (in pC/cm 2 ). 


Phase 


S.G. A E Gap 


T/ ac (A-AFM) 
T c (A-AFM) 
Taac (A-AFM) 
R aac (A-AFM) 
p (A-AFM) 
SG (FM) 


Pc 

Pna 2i 
Cc 
Cc 

Pnma 

C2/c 


138 

102 

93 

26 

15 

0 


2.1 

1.8 

1.8 

2.7 

2.7 

1.7 


(42,42,95) 
(0,0,75) 
(43, 43, 73) 
(61,61,39) 
( 0 , 0 , 0 ) 
( 0 , 0 , 0 ) 


FIG. 1. (Color online.) Energies of the bulk BiMnOs phases 
found to be either saddle points of the energy surface ( T' c and 
T ac ) or local minima (the rest), computed using PBEsol+C/ 
(top) and HSE06 (bottom), for each of four magnetic order¬ 
ings. 


ues gave good agreement between the structural param¬ 
eters computed for the bulk ground state and the experi¬ 
mental results; however, the obtained band gap was close 
to zero, and reasonable variations of U and J could not 
get the band gap to open above 0.5 eV (the experimental 
gap of FM BiMn0 3 is 0.9 eV^ but it is expected to be 
larger for antiferromagnetic configurations). Because of 
this failure to open band gaps, in many cases when opti¬ 
mizing FM structures we would get a spureous metallic 
state; this is the reason why in Fig. [T] (top) many of the 
FM states are absent. Another spureous result obtained 
when using PBEsol+G is that the experimental ground 
state of BiMn0 3 is predicted to have higher energy than 
the p phases. 

We then used the more computationally demanding 
HSE06 hybrid method to re-optimize the structures re¬ 
ported in Fig. [T] (top), and similar FM structures; for all 
structures but some of the GS ones it is possible to do 
this with 20-atom unit cells. All phases were found to 
be robust insulators, and the C2/c (GS) structure was 
accurately predicted as the one with the lowest energy. 
Otherwise, the energy differences between most phases 
are similar to those obtained using PBEsol+LL This is in 
line with the previous finding for BiFe0 3 in which we re¬ 
ported that energy differences between similar structures 
are independent of the exchange-correlation functional 
used due to error cancelation^ In particular, PBEsol+t/ 
and HSE06 predict very similar energy differences be¬ 
tween different magnetic arrangements of the same struc¬ 
ture; this is why we have not recomputed other magnetic 
orderings of the GS phase (some of them require to dou¬ 
ble the unit cell, making these calculations computation¬ 
ally expensive). For each of the geometries that are local 
minima of the energy, Tables Q] and El list in detail the 
properties of the structure having the favored magnetic 
ordering; in particular, the lattice parameters predicted 


for the C2/c phase are within 0.4% of the experimen¬ 
tal values^ The T phases have similar characteristics 
to those found for BiFe0 3 (an analysis of most of their 
atomic displacement patterns can be found in Ref. l25l) ; 
the most stable one has also a Cc space group, a large 
c/a ratio (1.17), and a large macroscopic polarization 
(94 pC/cm 2 ). Magnetic arrangements other than those 
shown in these Tables display similar structural configu¬ 
rations, differences in their lattice parameters being al¬ 
ways below 1%. The band gaps are always smallest for 
the FM ordering, and largest for the G-AFM orderings, 
the difference for a given structure going up to 1.5 eV. 


B. Structures of BiMn0 3 Epitaxial Films 

The bulk phases described here can be made stable 
by growing the material as a coherent epitaxial thin 
film. Like in BiFe0 3 j2£ the supertetragonal phases are 
expected to be favored at compressive strains; at ten¬ 
sile strains, the large in-plane lattice parameters of the 
p and Raac phases hint that these might be more sta¬ 
ble than the GS phase. To check these hypothesis we 
have done structural optimizations for the four lowest- 
energy phases, where the epitaxial effect is simulated by 
constraining the in-plane lattice vectors to be equal in 
length and to form a 90° angle. Figure [2] shows how the 
properties of these films change with the in-plane lattice 
constant. At lattice constants around 3.90 A the epitax¬ 
ial distortion on the GS film breaks the C2/c symmetry, 
but it is otherwise very small, so the corresponding film 
is still the favored one. As we compress the film, the T aac 
phase becomes competitive, and below around 3.75 A it 
is expected to be the stable state of the material. For 
tensile strains at in-plane lattice parameters above 4 A 
we expect the i? QQC film to be the ground state, displaying 
a large polarization with components both in plane and 
out of plane. The p phase is, by far, the one that has to 
be most distorted to fit the square symmetry of the sub¬ 
strate, and this renders it energetically not competitive 
with the other ones. The films are always insulating and 
display magnetic moments of around 4 pb localized in the 
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TABLE II. Lattice parameters and Wickoff positions of the 
bulk BiMnC >3 phases found to be a minimum of the energy in 
this work. 


Phase 



Structure 


Taac 

a = 

4.597 A, b = 5.230 A, c = 

= 5.256 A 



a = 

= 90°, /3 = 

91.3°, 7 = 

CD 

O 

O 


Mn 

2 a 

0.5667 

0.2552 

-0.0557 


Bi 

2 a 

0 

0.7799 

0 


O 

2 a 

0.1270 

0.8051 

0.4122 


O 

2 a 

0.6226 

0.5461 

0.1486 


O 

2 a 

0.6709 

-0.0276 

0.7141 

T a 

a = 

5.315 A, b = 5.214 A, c = 

= 8.764 A 



a - 

= 90°, P = 

II 

CD 

O 

o 

II 

90° 


Mn 

4a 

0 

-0.0046 

0.2788 


Bi 

4a 

0.5556 

- 0.0111 

0 


O 

4a 

-0.0303 

-0.0750 

0.0429 


O 

4a 

0.8105 

0.3081 

0.2885 


O 

4a 

0.7120 

0.7921 

0.3282 

Taac 

a = 

10.235 A, b = 

5.233 A, c 

= 5.310 A 



a = 

CD 

O 

O 

II 

121 . 2 °, 7 = 

= 90° 


Mn 

4a 

0.2210 

0.2546 

0.1668 


Bi 

4a 

1/2 

0.2602 

0 


O 

4a 

0.4572 

0.3218 

0.3688 


O 

4a 

0.7124 

0.4440 

0.3471 


O 

4a 

0.1703 

0.4634 

0.3953 

Raac 

a = 

9.266 A, b = 5.720 A, c = 

= 5.675 A 



a = 

II 

< 0 . 

o 

O 

125.7°, 7 = 

= 90° 


Mn 

4a 

0.5233 

0.2458 

0.3328 


Bi 

4a 

0.2422 

0.2422 

0.4844 


O 

4a 

0.2744 

0.3161 

0.1236 


O 

4a 

-0.0039 

0.4427 

0.0890 


O 

4a 

0.5754 

0.4725 

0.1539 

p 

a = 

5.928 A, b = 7.440 A, c = 

= 5.372 A 



a - 

= 90°, P = 

II 

C^~ 

O 

o 

05 

II 

O 

O 

C5 


Mn 

4b 

0 

0 

1/2 


Bi 

4c 

0.07245 

1/4 

0.01010 


O 

4c 

-0.02765 

i 1/4 

0.59634 


O 

8 d 

0.82708 

0.45595 

0.20851 

GS 

a = 

9.544 A, b = 5.592 A, c = 

= 9.852 A 



a = 

II 

< 0 . 

o 

O 

O} 

110 . 8 °, 7 = 

= 90° 


Mn 

4e 

0 

0.2921 

1/4 


Mn 

4c 

1/4 

1/4 

0 


Bi 

8 f 

0.6353 

0.2270 

0.1215 


O 

8 f 

0.5977 

0.1733 

0.5804 


O 

8 f 

0.1455 

0.0718 

0.3737 


O 

8 f 

0.3518 

0.0472 

0.1652 


Mn ions. Their c/a ratios grow markedly as the strain 
becomes more compressive, and this translates partially 
into larger out-of-plane polarizations. Most of the Mn-0 
bonds stay at values of around 1.9 A, but we can also 
see in Fig. [2] much longer bonds; those arise due to the 
orbital ordering in the GS phase, and in the out-of-plane 
directions of the T phases. 



FIG. 2. (Color online.) Properties of BiMn 03 films as a 
function of the in-plane lattice parameter: (a) energy relative 
to the bulk ground state; (b) band gap; (c) magnetic moment 
of the Mn ions; (d) macroscopic polarization; (e) c/a ratio; 
and (f) Mn-O distances. The four structures mentioned in the 
legend of (c) have been considered in their favored magnetic 
ordering (FM for GS and A-AFM for the rest). 


C. Magnetic Properties 

The different magnetic orderings of a particular 
BiMnC >3 structure are within around 50 meV/f.u., both 
in bulk (Fig. [T]) and in films (Fig.[3ta), top panels). The 
A-AFM ordering is favored both for the T acc and for the 
Raac phases. These phases cannot accommodate the or¬ 
bital ordering of the bulk, so the ferromagnetic ordering 
is not the prefered one any more. 

To estimate the Neel temperature for the T aac and 
Raac phases we used a Heisenberg model with energy 
E = E 0 + 1/2 Jij S i ■ S j ; E 0 is a reference energy, and 

Jij is the exchange coupling constant between the spins 
localized at Mn ions i and j , given by S; and (taken 
as unit vectors). We restrict ourselves to first neighbour 
Mn-Mn interactions in plane (described by J a ) and out 
of plane (described by J c ). Fitting our first-principles re¬ 
sults to this simple model we obtain the values for those 
constants that are displayed in Fig. [3](a) (bottom pan¬ 
els); the energies given by the model are represented by 
lines in the top panel, and they show reasonable agree¬ 
ment with the first-principle data. The exchange cou¬ 
pling constants for BiMnC >3 are very different from those 
of BiFeC> 3 ,— reflecting the different natures of Fe +3 (with 
a d 5 electronic configuration of half-filled orbitals, leading 
to strong antiferromagnetic interactions according to the 
Goodenough-Kanamori rules of superexchang o 30 : 31 ) and 
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Mn +3 (with a d A configuration that includes an empty e g 
orbital, favoring in-plane ferromagnetism—the Mn-O- 
Mn angles in the T aac and R aac phases are between 150° 
and 160°). We then used a Monte Carlo method to solve 
our Heisenberg model in a periodically-repeated box with 
20 x 20 x 20 spins; the results obtained for the order 
parameter that describes the A-AFM alignment and for 
the magnetic susceptibility are shown in Fig. [31(b). Bulk 
BiMnC>3 orders magnetically below around 105 and 
the Neel temperatures of our simulated films are similarly 
low—around 80 K for the T aac film with lattice param¬ 
eter a- m = 3.70 A, and around 90 K for the R aac him 
with a i n = 4.08 A. The exchange coupling constants do 
not change much in the ranges of epitaxial strains where 
the T aac and R aac are expected to be stable, so the Neel 
temperature will be similar in these ranges. 


D. Optical Properties 

BiMnOa is a particular perovskite oxide not only in its 
ferromagnetic properties, but also in its small band gap— 
typically, band gaps are above 3 eV for these materials,— 
but they are around 2 eV or less for some of the phases 
of BiMnC>3 according to our calculations and to the few 
experiments available— In order to further investigate 
the optical properties of BiMnC>3, we used the indepen¬ 
dent particle approximation implemented in Vasr^ to 
compute the frequency-dependent dielectric matrix, and 
the related absorption coefficient. Since no experimen¬ 
tal data for comparison are available for BiMnC>3, we 
did a initial test on BiFeOs—we compared the absorp¬ 
tion coefficient of bulk BiFeC>3 computed following this 
methodology with the one measured by Chen et al. for a 
film with a very similar structure— Our BiFeOs results 
in Fig.[3Ka) show good agreement between theory and ex¬ 
periment, if we correct a shift associated to the overesti¬ 
mation of the band gap by the theoretical method. (This 
band gap for BiFeOs (3.4 eV) is the same reported earlier 
by Stroppa and Picozzi^ who performed a calculation 
similar to ours.) For the BiMnOs T aac film expected to 
be stable when grown on a substrate of around 3.70 A, we 
find an absorption spectra that matches the solar range 
better than what is typical in other perovskite oxides, 
as shown in Fig. SKb). We found a band gap of around 
2 eV, but this could be even smaller if the HSE06 hybrid 
method is again overestimating it. This makes superte- 
tragonal BiMnC>3 films interesting in the framework of 
materials for photovoltaic devices where light absorption 
can be coupled to other functional properties. 




3.69 3.72 3.75 3.78 
In-plane parameter (A) 


3.99 4.02 4.05 4.08 
In-plane parameter (A) 



(b) 


FIG. 3. (Color online.) (a) Top panels: energy of different 
magnetic arrangements with respect to the most stable one 
(A-AFM), for the T aa c (left) and R aa c (right) phases. Bottom 
panels: exchange constants J from fitting those energies to a 
simple Heisenberg model, for the T aac (left) and R aa c (right) 
phases, (b) Top panels: A-AFM order parameter as a func¬ 
tion of temperature for a T aac film of in-plane lattice constant 
of 3.70 A (left) and for a R aac film of in-plane lattice constant 
of 4.08 A (right). Bottom panels: total magnetic susceptibil¬ 
ities for the same films (filled symbols) and their partial con¬ 
tributions, parallel (||) and perpendicular (_L) to the direction 
of the A-AFM order parameter (empty symbols). 


IV. CONCLUSIONS 

In this article we have explored the epitaxial phase di¬ 
agram of BiMnOs films with the help of first-principles 
calculations. While bulk BiMnOs is a paraelectric, we 
predict that it will transform to a supertetragonal phase 


with a polarization of around 100 pC/cm 2 when grown 
on substrates that compress its in-plane lattice constant 
to about 3.75 A (for example, YAIO3 or LaSrAlO^ii). 
This polarization will point mostly out of plane, but it is 
also possible for BiMnOs films to develop a polarization 
of similar size laying mostly in plane by growing it at 
tensile strains on top of substrates that expand its lat- 
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BFO ( R3c phase) 



BiMn0 3 (T aac phase) 



(b) 

FIG. 4. (Color online.) Imaginary component of the average 
of the diagonal elements of the dielectric matrix (left) and 
absorption coefficient (right) for (a) bulk BiFeC> 3 , and (b) 
a Taac film of BiMnOs (a; n = 3.70 A). The experimental 
values for BiFeOs were taken from Ref. |3J. The shaded area 
corresponds to the solar spectrum. 


tice constant beyond 4 A (such as BaTiC >3 or PZT— ). 
Our findings might explain the experimental reports of 


ferroelectricity in BiMnOs films grown on SrTiOa and 
LaAlOs with relatively large remnant polarization^— 
even if SrTiOs and LaAlOs have a larger lattice con¬ 
stant than the one needed for stabilizing the supertetrag- 
onal phase according to our calculations, these films show 
in experiments nonuniform strain distributions^ which 
could be the signature of a co-existence of different (po¬ 
lar and non-polar) phases of BiMnOs films. 

The supertetragonal phase that we have found is very 
similar to that of BiFeOs films grown on LaAlOj^ 
(ai n ~ 3.79 A) and of BiCoOs, where it is the ground 
state even in bulk (ai n ~ 3.73 JUft). This kind of phases 
with giant polarization might be ubiquitous in transition- 
metal perovskite oxides containing bismuth, whose lone 
electron pair can be easily accommodated in this kind of 
structures^ note that previous calculations for BiScOs^ 
provide a hint that this might be yet another material 
displaying the same behaviour. Since BiMnOs is the 
only strong ferromagnet among the insulating transition- 
metal perovskite oxides, our results are also relevant in 
the context of discovery of new multiferroics; we have 
found that A-AFM ordering is energetically favored, but 
the FM ordering is competitive (especially for the R aa c 
phase). It might thus be possible to engineer ferroelec¬ 
tric ferromagnets by combining these new film phases of 
BiMnOs with other transition-metals oxides in superlat¬ 
tices or solid solutions. Finally, in addition to this func¬ 
tional properties, some of these phases of BiMnOs show 
band gaps that are smaller than those found typically in 
perovskite oxides, which makes this material interesting 
also from the point of view of photovoltaic applications. 
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